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Electromagnetic production of pions and quark dynamical mass in Minkowski space 
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^Department of Theoretical Physics, Institute of Nuclear Physics Rez near Prague, CAS, Czech Republic 

Mechanism for generation of vector meson resonances is studied in the framework of QCD Dyson- 
Schwinger equations, which are dehned and actually solved in Minkowski space. It is suggested 
that the timelike pion form factor is generated by the interference of the background quark loops 
together with the resonant structure predominantly created in the photon-quark-antiquark vertex. 

It is suggested that QCD Green’s functions involving quark fields are oscillating for timelike argu¬ 
ments, which makes the interference effect among various QCD Green’s functions quite strong and 
important for the correct description of production processes. A further peculiarities of dynamical 
chiral symmetry breaking and confinement phenomena as viewed in Minkowski space are discussed. 

PACS numbers: 11.10.St, ll.15.Tk 


I. INTRODUCTION 

Probing the charged quarks inside of the hadrons by photons can provide nontrivial information about quark four 
momentum distribution in hadrons. Resulting charged pion form factor as well as of the transition pion form factor 
represent rather precise data on the electromagnetic structure of the light meson. In the spacelike region the available 
data [IHi are represented by a smooth decreasing line, however a certain discrepancy with asymptotic predictions 
Qbased on the perturbative QCD is exhibited. These day available experimental range ~ lOOGeU^ moves the 
perturbative QCD towards deep spacelike scale. 

Contrary to the spacelike behavior, the pion form factor iG in the timelike region represents the amplitudes for the 
production of 7r+7r_ pair with appearance of the striking resonant structure, it has a peak at 770MeV- the rho meson 
with a small (~ 5%) admixture of isoscalar oj resonance. Other resonant structure, e.g. the deep dip at 1.5GeV has 
been found by using Initial Photon State method at BABAR 2012 experiment where also higher energy peaks 
are observed. Recall here, that the pion form factor can be considered as a boundary value of analytical function, 
which has a cut on the timelike axis of square of momenta starting in the branch point Sth = Qth ~ 4mJ- the 
production threshold. On the other side, it is well known that the marriage of PT QCD and dispersion relation 
method is not possible without a conflict with requirement of the reflection positivity. In other words SP’tt cannot 
be positive function as sometimes assumed and the known serious attempts to evaluate dispersion relation should 
explicitly exhibit aforementioned violation of reflection positivity Q- It is notable, that also precise form factor fit 
made in Q does not respect a naive perturbative unitarity. Actually, the complex couplings are used for excited radial 
states of pion, which spoils the well known relation between the real and the imaginary parts of 

A reliable low energy description of QCD production processes, which is based on the QCD degrees of freedom 
was basically lacking and unknown. In Spectral Quark Model (SQM) [M3| a clever mechanism employed, if correct, 
is just the opposite; the cuts in the quark propagators conspire to build a pole and threshold in the form factor. 
Matching together hadronic and quark degrees of freedom is a typical feature of the models [m, which also try to 
explain behavior of TV at large Minkowski region of Q^. However as in the previous SQM, the dynamics is not 
obviously governed by QCD gap equations and supporting solution is waiting to be explored. 

The Dyson-Schwinger and Bethe-Salpeter equations (DSBSEs) are quantum equations of motions at the level of 
Green’s functions. This prominent framework of DSBSEs, which provides a traditional tool in the study hadronic 
physics, is up to a few exceptions implemented with the Euclidean metric due to the known numerical reasons. In QCD, 
using an Euclidean metric as a definite one, a various approximation of DSBSEs system provided a limited ground 
for hadronic physics, especially for determination of spectra of ground state mesons, their strong and electroweak 
decays and the spacelike (Euclidean) form factors |13 - [lq . In certain belief, similar has been attempted to extract 
information about excited mesons, however little is known about the implementation of the BSE for production 
process, where excited states are identified with broad resonances seen experimentally. One is obviously faced with 
obstacles of technical character: e.g. how to perform the analytical continuation of QCD Green’s function to the large 
timelike arguments, however also with more principal ambiguities: how to identify excited states with a snapshot of 
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resonances, especially when they are tricky hidden in the physical form factors. Roughly say, the pion electromagnetic 
form factor and associated PDG list of possible p and w vector resonance participants is the matter. For a practical 
deal with analytical continuation problem see |17H^ . for recent results for radial excitations based on solution of 
homogeneous BSE see for instance [2ll - l^ . 

In presented paper we continue the studyof QCD DSBSEs in Minkowski space, i.e. we are going to ignore 
the conventional wisdom and does not solve the DSBSEs in the Euclidean space, but instead of, the calculation is 
performed directly in Minkowski. To show that such, at least one, solution exists, is the main motivations of the 
author. Eor a good QCD DSBSEs model one should get light pseudogouldsteone boson of broken chiral symmetry- 
140 MeV heavy pion. Eurthermore, as very related consequence, one also should not get stable states with masses 
exceeding of the mass of several pions. In words, there should not be solution of homogeneous BSE for a heavy real 
mass for all p’s , scalars, radially excited pseudoscalars etc.. Instead of, all these broad resonances could be visible in 
a full color neutral vectors, scalar, pseudoscalar, axialvector, tensor... inhomogeneous BSE. Incorporating the ’’meson 
widths” in QCD DSBSE formalism means to consider nonperturbative solution of inhomogeneous BSE. We consider 
simplified solution of inhomogeneous vector BSE and use it for the calculation of pion form factor in Minkowski space. 

To achieve the convergence of the DSBSEs system, one has not only to ensure that the system is without uncontrolled 
singularities in all Minkowski space but also one needs to use a right choice of variables. Recall, a naive discretization 
of Minkowski space momentum integrals leads notoriously to not well defined numerics. Hence to solve the DSBSEs 
system, for which the method of the iterations is chosen, we take a maximum number of suited Lorentz invariant 
variable, leaving thus no space for badly convergent numerics. Similar method has been already successful in search 
of excited states in toy DSBSE with Minkowski metric [IJ, as well as in the Euclidean theory pH |. 


II. ELECTROMAGNETIC PION FORM FACTOR FOR TIMELIKE ARGUMENT 


The calculations of exclusive hadronic form factors, which are based on the quark and gluonic degrees of freedom 
were not available due to the complicated nonperturbative nature of low energy QCD. Let us mention the continuation 
of the perturbative QCD asymptotic behavior [26l - [2^ 


FTrit) —>■ 


647rV^ 

(II — 2/3n/) t Int/Agp^ 


( 2 . 1 ) 


derived for t = — —>■ oo to the experimentally available region does not reproduce the data for any 

Obviously, one of the main difficulties is the inclusion of resonances at the level of QCD CE’s. Observing the 
BABAR data resonances are spread over the all experimentally known range of s (s = >0). Another, 

quite important source of aforementioned fail has mathematical origin. In most popular nonperturbative approaches 
like lattice theory or DSBSEs here, the Euclidean metric is employed from very beginning, however the analytical 
continuation known from perturbation theory does not apply here. 

Actually all sophisticated non-perturbative methods in QCD are applied o nly in the Euclidean space. Thus for 
instance DSBSEs provides reasonable meson masse^ meson form factors [2^, as well as baryon electromagnetic 
form factor has been studied in spacelike domain Easy and simple models suffice to describe smooth form 

factor for < 0 since there are little traces of complicated resonant structure known from the behavior at > 0. 

The shape line of the p meson peak- is the most striking property of the pion form factor around s = (O.SGeU)^. 
Hence a various effective Vector Meson Dominance (VMD) models takes the mass rup = ITaMEV as a granted input 
parameter of a phenomenological effective Lagrangean. Eurther inclusion of heavier resonances is less certain game. 
Let us recall here the PDC 2015 list for p(I450) where the Breit-Wigner (BW) masses are spread from the value 
1265 ± 75MeV for e+e” —>■ tt+tt” [s^, passing one of the BABAR 2012 data fit, where this first excited p has a BW 
mass 1493 ± 15MeV for e+e“ —>■ 7r+7r“7 in Q and finishes with the highest PDC accepted value 1582 ± 17 ± 25MeV 

for e+e->■ ttottoJ reaction. Very similar, process dependent list, can be found in PDC 2015 for /9(1700) (the second 

p excitation this year), recalling highest BABAR fit, which is 1861 ± 17 in Q, being fairly apart from r lepton decay 
p excitation 1728 ± 17 ± 89 and being thus lOOMeV higher then the older fits based on e+e“ —>■ tt+tt" data 
[3l,[3l,|33. In fact, apart of the ground state, the BW mass values for 1 , isospin 1, light flawoured states listed in 

PDC are almost continuously spread over the all experimentally available energy region, some of them belong to the 
one state, some of them to another one. 

QCD dual VMD models are effective quantum field theory which deal with a strongly coupled rho meson fields. 
The classical mass of given vector is identified with some central value of experimentally determined BW parameter. 
The excitation of spin one fields become a wide resonances due to the quantum effect- in practise , after a single loop 
is calculated. This large widths of all p’s makes their labeling by radial quantum number quite puzzling. We will show 
that this is a large interference effect between various QCD CF’s which makes the situation even more cumbersome. 
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In the expressions for form factor, the resonances, which originates as a poles in the quark-photon vertex are further 
weighted by the function which are oscillating in the timelike region. While it does not offer an easier and intuitive 
picture, the resulting form factor can be calculated in practice. It is a matter of the fact that the first numerical 
results suggest that there is no one to one correspondence between PDG (VMD modeled) list of resonances and the 
one we need to use in Green’s functions description in our DSBSEs model. Furthermore, not each DSBSEs resonance 
need to have its counterpartner in phenomenological VMD (and even in PDG). It will be exemplified explicitly for 
the case of . 

The evaluation of the pion form factor is a typical quantum field theory problem and how to calculate such transition 
in BSE approach is generally known [s^. As an additional constraint we should respect (axial) Ward identities- the 
consequence of (chiral) gauge invariance expressed in the terms of Green’s functions. At least one such an expansion is 
known [mu and here we will consider only the first term, which defines the so called (dressed) Relativistic Impulse 
Approximation (RIA) 

For the case of electro-production of pions the RIA reads: 


G'^(P+,P_) = eF^{q^)iP+- P.y 
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where in the first line the photon with momentum Q = P-^+P- couples to the up-quark with electric charge 2/3e while 
the second line describes triangle diagram where the photon couples to the down-quark. The function P bs stands for 
negatively charged outgoing pion,i.e. for BS vertex function with the total momentum P_ = —P + Q12 and relative 
momentum given by the difference of momenta fci and k^. The matrix P^j^/ is the quark-photon vertex determined 
by its own inhomogeneous BSE (shown bellow (12.101) '). Analogously we have P+ = P + Q/2, the quark propagator S 
connecting positive and negative pion is chosen as k^ = k — P, the propagators connecting electromagnetic vertex are 
ki ^2 = kF Q/2. In the isospin limit the propagators of u and d as well as the quark-photon vertices of u and d quarks 
are identical. By applying charge conjugation one can show that the second line, up to the different prefactor, which 
turns to be -|-l/3e, is equal to the first one. 

All the functions entering expression (12.21) should be self-consistently obtained from the solution of QGD DSBSEs 
system. Neglecting electroweak interaction, the quark propagator S can be readily written 


S{p) = [A{p) B{p)] 1 


1 i> + M{p) 

Ap"^ — M^{p) ’ 


(2.3) 


where the function M is usually called dynamical quark mass and the function A is the inverse of quark renormalization 
function Z. In the next section we present the details of QGD DSBSEs model, which not only provides the correct 
pion spectrum but which, after substituting the results into Eq. (1^ provides correct electromagnetic pion form 
factor. 


A. Quark gap equation and pion BSE 

In the following the Minkowski metric is employed. The scalar product of two four-vectors is defined as p.q = 
with metric = diag{l, —1, —1, —1, —1), thus p^ < (>)0 for spacelike (timelike) p. 

The propagator satisfies the DSE (the gap equation), which reads in momentum space 

S~^{p) = Z2 i)- ZiUigip) - ZiE{p) 

S(p) = i j ^^^G^^(p-fc)y7^S'(fc)r“(p,g), (2.4) 

where m{p) is the renormalized quark mass at the scale p, r“(p, q) is the quark-gluon vertex satisfying its own SDE 
and G^j/ is gluon propagator. 

The BSE vertex function r^^(p, P) as well as the BSE wave function x(fc,P) = S{k+)T{k, P)S{k-) is solution of 
bound state BSE: 


f d^k 

J ( 2^)4 


K{k,p, P)x{k,P)i^ 


TiJp,P) = i 


(2.5) 
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where K{k^p,P) is the renormalized quark-antiquark interaction kernel and where the total momentum satisfies 
. Strictly speaking, the homogeneous BSE (12.51) represents narrow mass approximation at the resonance 
mass for P^ = , where is the mass of the pion (for n = 1) or of the arbitrary excitation 7r(1300), 7r(1800),... . 

The pion is well established pseudo-Goldstone boson of broken chiral symmetry in QCD. Independently on the 
metric, in order to get a correct properties of the pion a symmetry preserving truncation of infinite DSEs system is 
required. The simplest Ward-identity preserving truncation is the so called rainbow-ladder approximation (RLA), 
which we are going to use in our Minkowski space model as well. 

The RLA enables us to write down the effective charge a: 

-p)^ {k - p)a{k - p) (2.6) 

defined as a simplified product of the quark-gluon vertex and the gluon propagator. We use the same method of the 
solution as in the previous paper , however here we have included the renormgroup improved UV part as well. 
The kernel of BSE and DSE is defined as following 


K(x) 


K,f{x) 


E{x) 

G-\x) 


P-UV (^) 


x{x) 
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d r 


-0(-*)— exp(-£(-cr))i/2in(G-i(x)) 
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In [e -I- (x/a)^/^] 
(x — a^) -I- 

Inin 
dx 


H 2N 


^'■QCD 


(2.7) 


where x = {k — pY and the capital letters G stand for the effective infrared and ultraviolet coupling strengths. The 
later we fix through one loop perturbation theory from which one gets Cuv = ( 33 - 2 ^j.) 27 i- ~ 0.0714. The scale a 
of dimension [massY is a parameter of the model and is set up uniquely after the identification of the pion mass. 
It defines log modulated energy period of the kernel oscillations at the timelike regime and exponential 

weakening for a large spacelike momenta. 

Well known intrinsic property of the dynamical mass function M in Eq. (12.31) is its fast increase in the infrared, 
where the it reaches some typical constituent quark mass value M(0) ~ 2 to^. Here however, the mass function M is 
a complex object which phase appears to be sensitive to the ratio of infrared and ultraviolet coupling strength. Let 
us note, that the physical pole mass defined as S~^{m) = 0 arising in non-confining theory like QED, does not exists 
in QCD due to the confinement. Here it is assured by the solution, M becomes complex and oscillating function with 
imaginary part never small at vicinity of such p where p = ReM{p). To achieve numerical convergence, the correct 
value of I3qcd is reached step by step by replacing 

Cuv kCuv ) 


where k is a small auxiliary number introduced at the beginning of the iterations and finishes with the value k = 1 
without lost of numerical stability during the iteration process. 

More interesting is the behavior of the function M itself, for small k it has oscillating behavior due to the oscillating 
kernel. When Cuv strength is pronounced the frequency of oscillations is lowered in a special way. For small k 
the function M is fairly oscillating at the both - spacelike and timelike domain of the Minkowski momenta (see 
FiglTJ increasing k the oscillations turns to fluctuations at the spacelike domain of momenta where it finnaly leads 
to occurrence of several extrema in QCD limit k ~ 1. If the ratio Kuv/Kinf was free the details of such behavior 
would be dependent as well. However, this ratio is fixed by the property of pion (by mY and /tt values) and it seems 
to be rigid property of our model, that when we finally add the complete UV part in Eq. (I131)( i.e. when K = 1), the 
solution remains relatively slowly oscillating at the timelike region, while this oscillation completely disappear in the 
spacelike region of momenta. Assuming small freedom given by the choices of Kinf and a, the open possibility is the 
behavior of 5R(S)M with n = 1 ± 1 crosses of zero at negative axis of momenta p^, for actual solution see Fig. [5J 

Quenched, RLA offers the solution Z{p) = 1 in Landau gauge, while the other solutions in the Euclidean space 
typically give Z{0) ~ 1/2 when Z{oo) = 1 in Landau gauge. We have found the renormalization function Z(j)) should 













M/a 




FIG. 1: Solution of QCD gap equation with Kinf only. Left panel shows the quark dynamical mass function as a function of 
Minkowski square of momentum The same function and the analytical ht (EU are displayed on the right panel against Inq^. 
bf chybicka to stou UV c“ 



FIG. 2: Solution of QCD gap equation with the full kernel K. Here, ReB is negative and ImBiQ in infrared. 


not be neglected in our model, since it significantly controls the infrared value M{p 0). Its introduction costs the 
numerical stability at some parameters space, and we make a real valued approximative analytical fit of the following 
form 


Z{P^) 


l + Czlni/^(e + |pVA^| 


( 2 . 8 ) 


which is valid for the complete kernel case {k = 1) and where reasonably numerical fit gives the constant Cz = 2. 

The gap equation has been solved after the renormalization as well as it has been solved in its un-renormalized 
form. Using a Pauli-Vilars regulator and large spacelike renormalization scale we got basically the same results as 
for the un-renormalized case. Hence we freely take for renormalization constants Z\ = Zi = Z/^ = \ vt\ thorough the 
paper. To complete the model, we take the current quark mass rUc = m{oo) = 0.01a and Cif = 3/100 and UV par 
is fully taken into account. For completeness we should mention here, that there exists more then one solution of the 
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DSBSEs system and in many cases we found evidence for at least another nontrivial solution. To get rid of unphysical 
solutions, we have calculated BSE as well as the observable F^^F* from both of two numerically identified solutions, 
skipping here completely the discussion about the one, which does not lead to satisfactory description of the pion 
form factor. 

In order to get the meson spectrum, one first solves the gap equation (El, whose solution is required to complete 
BSE kernel, and then solve the BSE for pion and its radial excitation. Here we consider the dominant component of 
the BSE vertex function, i.e. 


r(/c,p).„ =75ri5,„(fc,p) 


(2.9) 


which is known as good approximation already from early studies of DSBSEs [40| . 

The quark propagators need to be evaluated for various arguments p± as it appears in the BSE kernel. In order to 
eliminate one of the sources of potential numerical instability, we avoid the interpolation of numerical data and we 
make an analytical fit of the DSE solution. The fits, which have been actually used are presented in the Appendix A. 
A further important facts about Minkowski numerics are discussed in the Section III. 

As mentioned the UV part has a large influence on the solution and in fact its presence also change the scale 
parameter a, which appears in Kinf. Actually, for the pion mass we get Mtt ~ 0.7a for /c = 0, while we obtain 
Mjr ~ 0.2a for /c = 1. Contrary to the similar Minkowski model 2^ we do not get the solutions for excited states 
here. By looking the eigenvalue of BSE, we see only the evidence for the first excited state, which is located at 
M = 1300 ± \00MeV. We argue, instead of homogeneous one, the inhomogeneous BSE must be solved to see a broad 
resonances like 7r(1300) and 7r(1800). 


B. Electromagnetic vertex 

The information about bound states is encoded in the solution of the following inhomogeneous BSE for the quark- 
photon vertex: 


rUb, P) = P)s{k-)K{k,p, p) . (2.10) 

The resonances correspond with the complex poles in Tem in total momenta P^. As there is no stable narrow 
bound state in this vector channel the resonances do not decouple for a physical, i.e. real valued momenta. The 
function P) appears linearly in (12.21) and thus vectors resonances represent most important contribution to 

the pion form factor. The spin 1 axial vectors and higher spin tensor resonances contribute to P^- as well, however 
corresponding diagrams are suppressed by agEO since at least two photons should be attached to the loop made out 
of the quark lines. They do not contribute in RIA and we will neglect their contributions here. We also ignore puj 
mixing. 

Summation over the color is implicit in Eq. (12.101) and the function K{k,p,P) is the four-quarks scattering kernel, 
which should be taken in the form defined above (12.61) in order to maintain char ge c onservation. The solution of (I2.10|) 
consists from the twelve tensors accompanied with twelve scalar form factors [4ll |. The four of them are driven by 
Ward-Takahashi identity and can be called longitudinal part of the vertex T, the eight remainders are transverse and 
involves pre-formation a rho meson peak. After the suitable projection, the BSE turns to coupled set of 12 equations 
( in practice it leaves us with 8 equations as the four form factors are fully determined by QED WTI). 

Minkowski space solution of (12.101) remains basically unknown. The author was able to solve the system of equations, 
where in addition to known longitudinal terms, just the single following transverse component 

IrFsiQ^p) = ( 7 ^ - ^^)F3(Q,p) 

was considered (we keep notation [H). It has provided a rather complicated behavior of the function F 3 {Q,p) which is 
spanned above three Lorentz scalars (convenient choice can be associated with fermion legs momenta = {q±Q(2)‘^ 
and the photon square of 4-momentum Q^). We found that in the model defined by Eqs. (12.lOp . (1^^ . (12.611 . (12.71) the 
function P 3 is oscillating whenever one of the arguments is timelike. Unhappily, the odd attempts to include other 
Lorentz structure shows the equation system suffer by numerical instabilities and is not yet in the stage which can 
provide benchmark for precise evaluation of the function F^^. Due to this unpleasant fact, we will use a mean value 
Anstaze for single component approximation of Eq. (12.101) . For this purpose we take 


rU(9±,0') = 7'^(i + r.(Q')) 


(2.11) 
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i.e. the averaging is formally taken over the arguments of the quark-antiquark square of momenta and where newly 
introduced function satisfies rr(0) = 0 in order to maintain correct pion charge normalization = 1. Obviously, 

the electromagnetic vertex does not satisfy WTI, as we have neglected relevant longitudinal terms. While important 
from theoretical reasoning, these terms gives a second order effect when compared with p contribution and we ignore 
them completely thorough this paper. 

Using the tree level diagram for represents the first simplest step. In this, let us call it -Background Relativistic 

Impulse Approximation- (BRIA) one ignores the resonances completely by taking Tr[Q'^) = 0, i.e. = 7 ^. 

Obviously BRIA is not satisfactory and for purpose of the next discussion, we review some details of the Fpi fit 
made by BABAR collaboration itself @ . The fit is based on Gounaris-Sakurai VMD model [l^ , where the pion form 
factor F^ is made out of the four p BW resonances with masses mp^{785, 1493,1861,2254) where each resonance is 
characterized by its own complex coupling strength. We refer their absolute values 

\ci\ = (1,0.158,0.068,0.0051) (2.12) 


for completeness here. 

As we will see the DSBSEs model we use, will naturally explain these quite distinct couplings as a consequence 
of the interference effect among quark propagators, the electromagnetic vertex (12.101) and pion vertices as well. The 
effect of the interference among GF’s is abrupt, and as mentioned, the Ansatz for BSE solution, which comply the 
data, differs from the BABAR ht. The ’’spectrum” for the vertex BSE, which works phenomenologically is quite close 
to the (relativistic) harmonic oscillator 


~ i 

instead of string inspired Regge trajectory Mf ~ F To mimic resonances, we use the following form factor 


F,(g2) 

nQ) 


n _ Q'^ 

^ ml 




(2.13) 


where the functional form of the single phase (j)EM is taken to be similar as in the quark-gluon vertex (lO) . Explicitly 
written the phase we use reads 


4>em{Q‘^) 


(p j a? 


In (e -I- q^/a?)) 


(2.14) 


where /e is the constant prefactor. In other words, we assume that both gauge vertices are oscillating in the timelike 
region however their functional forms slightly differ there. Contrary to pure QCD quark-gluon vertex, the color singlet 
quark-photon vertex includes sign changing enhancements which turns to form factor resonances at the end. The two 
real constants mp ^, 7 p^ are enough for each individual BSE resonance, no special prefactor (the analogue of (|2. 121) 1 will 
be needed to explain decreasing crossection at higher . The list of parameters fit is in Tab. I. It is not completely 
surprising fact, that the list of resonances based on the vector BSE fit can differ from the one based on VMD (and 
hence PDG as well). For instance, the fit used here does not prefer the 1.5GeU heavy resonance used tho generate 
the pion form factor dip in VMD and there are two lighter vertex resonances instead. This is a certain and notable 
phenomenological difference between a various VMD and DSBSEs model presented here. 


III. NUMERICAL SOLUTION 


It is well known empirical fact that a naive integration over the space-time components of the loop momenta does 
not lead to a convergent solution. This numerical badness is a property of many tested model and is common to the 
systems which have regular and ultraviolate finite kernels as well. Published DSBSEs studies of the pion form factor 
1^ and the spectator approximation ^ are challenging, but employed for spacelike (Euclidean) only. A 
novelty, presented here, is the integration of DSBSEs system, which is performed directly in the Minkowski space. 

Remind for clarity the procedure, as a first one needs to solve gap equation, then one have to look what are the 
outcomes for the pion by solving the pseudoscalar homogeneous BSE, after that, when one knows the ratio m-^ja 
one can calculate the pion form factor. This procedure is usually repeated several times as one needs to improve the 












parameters of phenomenological fit (I2.13L which in our case finishes with the values presented in the table. Here we 
describe all details, without which the numerical solution would be actually impossible. 

We use LRA BSDSEs kernel with known primitive function by our construction, which not only improves accuracy 
but significantly shrink the time of iteration procedure used in our numerical search After an easy substitution, 
the angular integration over 3d subspace was done analytically, which leave us with two numerical integrations in 
DSE as well as in the BSE case. In this way we have avoided the use of angular expansion, which while working in 
the smooth numerics of Euclidean space, would turn to inconvenient method in the case of Minkowski space. 

Further , in order to integrate vector BSE (not shown here), and especially to integrate pion form factor (12.21) we 
always use Lorentz invariants for the numerical integration. Their concrete choice is quite arbitrary as follows from 
translation invariance. We have found that taking one, say the loop momentum is actually efficient for our system 
of gap/pion equations. However to find convergent integrator for Eq. (12.21) requires more. In this case we have found 
advantageous to identify integration variables with the arguments of the quark propagators for this purpose. 

First let us choose the c.m.s. frame defined by the photon momentum = (Q, 0,0,0), then we can take x = 
{Q/2 — k)‘^ ; y = ((5/2 + fc)^, where the inverse substitution gives — k —= •^/((fco —<5/2)^ —a: and ko = {y — x)l2Q — Ql2. 
In this way one directly gets for the Lorentz invariant measure the following expression 






(3.1) 


where D = — Q/2)^ . Outgoing pions are physical, i.e. their four momenta satisfy (p ± <5/2)^ = , which 

implies that relative momentum is purely spacelike. In our notation = (0,0, 0,pz) and for instance the product p.k 
is given as the following 

k.p = Q^/4 — cosd . (3.2) 

Actually, the formal mathematical operation expressed by the Eq. (13. ip defines new measure, where only the right 
hand side works numerically after the integrals are discretized to the sums (for which we use the Simpson method with 
Gaussian distribution of mesh points). While not needed here, let us note for completeness, that the Eq. (I2.10p has to 
be solved off-shell, thus an analogue pair of x, y internal variable, say: x’ = {p—Qf2)^, y’ = {p+Q/2)^ external variables 
are required for the introduction of the time component of external relative four-momenta Po = (y' — x')l2Q. Without 
aforementioned substitution one does not get convergent integrals in the equation for form factor and numerical noise 
with unwanted oscillations would prevent the stability of the iteration process employed (the all badness of using 
i.h.s. m can be seen by evaluating the form factor (j2.2l) . where obviously no iteration is needed). The reason why 
a naive integration over the Cartesian coordinates does not work, while the integration with Lorentz invariants works 
numerically remains unknown mystery, at least to the author. All relevant codes are available at the author’s web 
page (??). 

For the same reasons as in the case of the BSE solution we use the analytical fit for the all functions entering the 
integrand of For BS vertex we use chiral approximation r(p, P —>• 0) = M{p)/f^, improved due to the small 

nonzero current mass in the following way: r 7 r(p, = to^) = 
model. 

After evaluating the traces in (12.21) . the numerical integration is performed as explained above. Resulting form 
factor is shown in Fig. [3] for several values of fq and fe (see Eq. (12.141) 1. The approximative equality fq ~ fe is 
affirmating similarity of oscillating behavior in the both GF’s : in the quark-photon and in the quark-gluon vertex. 

We have found F does not posses spectral representation with 4 to^ threshold value, instead of, it is complex function 
everywhere (up to = 0 where the normalization condition P(0) = 1 ensures reality). Actually , after the correct 
normalization, the imaginary and real parts of F are not related via dispersion relation. Obviously, F is not direct 
observable in quantum field theory, but the scalar product FF* should be. At this pragmatic level, ignoring the 
question of analytical unitarity here, the same physics will be certainly achieved when 

FoSE/BSEy) = Fan{<f)F^''^‘^ ^ , (3.3) 

where is some running phase showing an allowed deviation from an ideal function Pq„, assuming Fan is the 

analytical function with the usual cut, and it exists. 




. where we take L = IGeV for our 


IV. SUMMARY, CONCLUSION AND PROSPECTS 

In this paper we have overcome usual numerical obstacles and solve the DSBSEs system directly in Minkowski 
space. To reach numerical stability, we use maximal number of Lorentz invariants as integration variables which turns 
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FIG. 3: Product F-^F^ calculated for several values fq,fe, which pair of parameters are used to distinguish the lines. BABAR 
data (crosses) and BRIA are also shown. 


usually unstable Minkowski numerics to the stable one. The numerics is at least enlightening and quite promising for 
future applications. The pion form factor has been calculated in RIA with phenomenological vector Bethe-Salpeter 
vertex contribution. It shows the importance of the interference effect between Bethe-Salpeter p mesons peaks and 
the background quark loop. The physics of resonances is quantum field phenomena where the interference effects do 
not play negligible role and can explains the complex couplings used in VMD models. However, the setup and shapes 
of resonances in DSBSEs approach is still based on phenomenology and a solution of more complete vector BSE is 
required in the future study. Presented numerical calculations shows oscillation as well as it suggest the complexity. 

Apart the BABAR fit, we did not list here a huge list of references devoted to pion form factors, however for a 
recent treatment of old-fashionable VMD model within the use the dispersion relation technique see 0,0,141, for 
the use of ChPT see for instance [i^, 0| and Regge approach used in the pion form factor calculations can be found 
in [isl - f^ . recall here also Resonance ChPT used in hadronic t decay Monte Carlo simulations [11. Most of 

the approaches do not contradict with analyticity, while the disagreement of QCD with spectrality is known for years. 
However here, if one is willing to accept the scenario proposed, the analytical structure of QCD GF’s requires further 
attention. As QCD Green’s functions here, are not an analytical continuations of their Euclidean counter-partners 
it leads us to the result for hadronic form factor, which does not comply with the dispersion relation as well as with 
Perturbative Theory Integral Representation [54l| . As a consequence, the pion form factor does not have spectral 
representation, and it suffers by the presence of unobservable ’’confining” phase. 

Very interesting would be a study of heavier mesons. Especially, for bottomonia we expect the aforementioned 
interference phenomena should be less pronounced due to the more ’’particle” nature of heavy b quarks. However even 
there, the so called hindered transitions [ssl - fs^ are very good experimental candidates, where overlaps of GF’s with 
a large interference effect can be expected. Various XYZ candidates [H^, reminding the famous state Ar(3872) [s^ 
observed by many experimental collaborations , shows up as enhancement appearing in more then one spin channel 
MM- Obviously, without necessity of fit to a simple quark-antiquark bound state picture, such enhancements or 
experimental bumps can simply originate in the interference effects of QCD Green’s functions. 


V. APPENDIX A 

The following fits of quark DSE solutions have been found useful for evaluation of electromagnetic form factors. 
They ht the solutions for different ks, hence different ratios m.n./a, where a is the interaction kernel parameter and 
niTr is the pion mass. 

For the case with no no mixing in DSBSEs kernels ((k = 0)), one gets mT^ja ~ 0.7 and the following fit has been 
found for the mass function: 
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ReM (x) = — 


S{x) 

0.1 + X 


cos uj(x) + mcZ(x) 


ImMix) =-sina;(a;) 

^ ^ 0.1 +a; ^ ^ 

L{x) = In [e + \/x/l-8 + a;/1.4 + (a;/30)^] 

'/x 


S(x) = aL(x)Z(x) ; uj{x) = fq 


L{x) 


(5.1) 


for spacelike momenta < 0, where we use shorthand notation x = —q^ja^ and the fit : 


ReM(a;) 

ImM (x) 
T{x) 


2T{x) 
1 + a; 


cos u}{x) — T{x) sin u}{x) 


QiT (x) 

— sma;(a:) — T{x) cosw(a;) 
1.7aIn^'^ (e + x ^), 


(5.2) 


which is valid for timelike momenta, where for now x = q^ ja?. For an achieved accuracy of the fit see Fig. [T] The 
numerical fit gives fq=8 for the kernel with((tt; = 0)) . 

For combined kernel with k = 1 one gets m.,^1 a ~ 0.2 and the spacelike part must be reconsidered. The following 
fit has been used for /c ~ 1: 


ReM{x) = ^ 

C2a + X 

lmM{x) =-— (5.3) 

^ ^ C2a + a:: ^ ^ 


for spacelike momenta q^ < 0, where typically we took C 2 — 0.01 — 1 and Ci ~ 0.02 — 2. 

The fit for the timelike M remains the same as in even for combined kernel k ~ 1 /2 — 1, since we have found 
that variation of fq was sufficient for our purpose. For (k ~ 1) we found the quark mass frequency lies in the interval 
fq = (2,3) very independently from the value Cinf used in the kernel. Combinations of these two fits ( 15 . 11 ) and ( 15 . 3L 
assuming simultaneous readjustment of the parameter fq provides high quality fit for any value of k , providing thus 
the solution in the range mT^ja ~ 0.1 — 0.7. Let us stress, the pion form factor has been evaluated for the real QCD 
strength, i.e. k = 1. 
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i 

1 

2 

3 

4 

5 

6 

m 

775 

1000 

1240 

1850 

2350 

2950 

7 

220 

250 

340 

250 

250 

350 

msfit 

775 

- 

1 493 

1861 

2254 

- 

mpDG 

775 

- 

1450 

1700 

2150 

- 


TABLE I: First two lines provide two parametric fit characteristic for six resonances of the model, ms fit are the values of 
complex BW fit made by BABAR. The last line are masses of p resonances presented in PDG 2016 , related with 7r''"7r“ mode 
in PDG (there are also p(1570) and p(1900) states in PDG which are not measured in two pions production). Instead of one, 
our model have two resonances 




